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We report a combined experimental and theoretical simulation of multiply scattered light dynam¬ 
ics in an ultracold gas of ®®Rb atoms. Measurements of the spectral dependence of the time-decay of 
the scattered light intensity, following pulsed excitation with near resonance radiation, reveals that 
the decay for long times is nearly exponential, with a decay constant that is largely independent 
of detuning from resonance. Monte Carlo simulations of the multiple scattering process show that, 
for large detunings, near resonance scattering of Fourier components of the excitation pulse plays 
a significant role in the effect. This interpretation is supported by the observations, and successful 
modelling, of beating between Rayleigh scattered light at the excitation carrier frequency with the 
Fourier components of the excitation pulse that overlap significantly with the atomic resonance. 

PACS numbers: 


I. INTRODUCTION 

Multiple scattering of near-resonance radiation is well 
known to have significant effects on the density distribu¬ 
tion of trapped ultracold atoms [ij , and on the collective 
dynamics of the trapped ensemble itself 0. However, in 
spite of the enormous number of applications and stud¬ 
ies of atoms confined in a MOT, there have been rela¬ 
tively few direct studies of the dynamics of the multi¬ 
ply scattered light itself. The limited number of studies 
include measurements of radiation trapping Q of res¬ 
onance radiation in an ultracold gas of atomic Cs by 
Fioretti, et al. and studies of the long-time dynam¬ 
ics of light scattered in the vicinity of the atomic ®^Rb 
F = 3 —>■ F' = 4 resonance transition ; these mea¬ 
surements were also made in an ultracold gas, for which 
there is limited frequency redistribution of the scattered 
light. Recent measurements by Balik, et al. Q examined 
the influence of quantum interference on the dynamics 
of the atomic alignment produced in the vicinity of the 
F = 3 —?► F' = 4 hyperfine transition associated with the 
atomic ®®Rb D2 resonance transition. 

Multiple light scattering is also important in several 
specific areas of current research interest [sl-fl^. These 
areas are focussed on the role of interferences in multiple 
wave scattering in a wide variety of condensed, liquid, 
or gaseous samples. Such interferences are readily ob¬ 
servable in many condensed systems and under a broad 
range of circumstances. For warm atomic gases, the in¬ 
terferences are normally not observable because of fre¬ 
quency redistribution caused by the thermal motion of 
the constituent atoms or molecules. However, in ultra¬ 
cold atomic gases [TMl , where the thermal motion of 
the atoms is very small, the interferences become read- 
ily observable through the coherent backscattering effect 
[iTHl^ . Research to date has recently been reviewed by 
Havey, et al. and by Kupriyanov, et al. [2l|. In rel¬ 
atively low density gases, where the mean free path for 


light scattering is much greater than the mean separation 
between the atoms, the process can be considered to be 
due to a sequence of separate scattering and propagation 
events. All studies to date have been done in the so-called 
weak localization regime, defined by the loffe-Regel con¬ 
dition as occurring when kl >> 1. Here k is the wave 
vector of the light, while I is the mean free path for light 
propagation in the medium. The coherent backscatter¬ 
ing enhancement is due to distributions of chains of such 
sequences with a corresponding distribution of scatter¬ 
ing orders. Scattering along reciprocal paths then leads 
to the interferometric enhancement of light scattered in 
the backscattering direction; in this direction the relative 
geometrical phase of the reciprocal paths is preserved, 
irrespective of the location of the atoms comprising the 
chains. When the density of atoms is increased to the 
point where the mean free path for light scattering is 
on the order of the atomic separation, then recurrent 
scattering Q becomes important, and a transition from 
propagating light modes to localized ones becomes possi¬ 
ble. The boundary near which this is expected to become 
important is defined by the loffe-Regel condition as kl ^ 
1. Localization by disorder (in this case disorder due to 
the random distribution of atoms in the sample) is the 
optical analog of Anderson localization [2^ of electrons. 
One driving motivation for studies of diffusive light prop¬ 
agation is to understand the purely diffusive part of the 
light propagation, and how the contrast with a transition 
to localized modes might be observed. Light localization 
has been reported for condensed samples in the optical 
M and in the microwave regimes [ll|, but not in an 
atomic gas. 

A second area of considerable interest is coherent ma¬ 
nipulation of light propagation in ultracold atomic gases 
1^ . This is accomplished by using effects associated 
broadly with coherent population trapping [23 - [^ . and 
specifically under conditions associated with electromag- 
netically induced transparency (FIT) [T7I - [30I| . One re- 
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markable demonstration of this was the seminal experi¬ 
ment by Han, et al. 31|, in which reduction of the speed 
of a coherent beam of resonance radiation to 17 m/s in 
ultracold atomic sodium was observed. In this experi¬ 
ment, the group velocity of a weak probe beam was ma¬ 
nipulated by a second control beam, where the two light 
sources were tuned to a two-photon Raman resonance 
between Na ground state hyperfine levels. Other exper¬ 
iments have used EIT to demonstrate nonlinear optical 
effects inclu ding four-wave mixing [s^ . optical informa¬ 
tion storage [33|, remarkable recent experiments demon¬ 
strating generation and control of single photons 
[ 3 ^, and general studies of EIT phenomenology in lad¬ 
der and lambda type configurations [s^- All such stud¬ 
ies have been concerned with single or multiple coherent 
beams either applied or generated by nonlinear processes 
in the samples. However, as demonstrated by the co¬ 
herent backscattering effect, phase can also be preserved 
in multiple scattering of light in an ultracold atomic gas 
[ 33 • This leads to the possibility that the CBS effect, 
in the weak localization regime, can be modified by the 
presence of a second control field in one of the customary 
EIT configurations. Recent theoretical results have in 
fact demonstrated that profound effects on the CBS in¬ 
tensity profile, spectral variation, and the temporal pulse 
shape of multiply scattered light beams are greatly mod¬ 
ified by the presence of a control beam [s^. A physical 
way to envision these effects is to realize that the CBS 
effect reflects the optical transport properties of a weak 
probe beam, and the presence of an additional control 
field modifies these transport properties. This modifica¬ 
tion of the diffuse light propagation in an EIT configu¬ 
ration is a newly discovered effect, and much exploration 
of the experimental and theoretical landscape remains to 
be done. 


II. EXPERIMENTAL APPROACH 

The experimental details of our approach have been 
described elsewhere Q , and will be only briefly sketched 
here. The basic experimental scheme is portrayed in Fig. 
1, where it is shown that physical sample is an ultracold 
gas of atomic ®^Rb prepared in a magneto optical trap. 
The trap operates on the F = 3 —^ F' = 4 hyperfine tran¬ 
sition. The trap [l^ produces a nearly Gaussian cloud of 
^ 10® ultracold rubidium atoms at a temperature ^ 100 
/iK. The peak density is ~ 3 x 10^° cm“®. Fluorescence 
imaging is used to estimate the Gaussian radius of the 
sample to be rg ~ 1mm. The peak optical depth, as mea¬ 
sured by absorption measurements of light transmitted 
through the central portion of the atomic sample, is bg = 
8(1). For a Gaussian atom distribution in the trap, the 
weak-field optical depth, on resonance and through the 
center of the trap, is given by bo = V^nocrovo, where 
no is the peak trap density and ctq is the on-resonance 
cross-section. The single resonance scattering cross sec¬ 
tion cr varies with probe frequency, b = 6o[l + ( 2 A/ 7 )^]“^, 



FIG. 1: Schematic diagram of the experimental arrangement. 
Shown in the figure is the probe laser, magneto optic trap 
(MOT), linear polarizers (LP), and a photomultiplier tube 
(PMT) used for detection of the fluorescence signals. 


where A = wl — wg, and wl is the probe frequency, wg is 
the F = 3 —>■ F' = 4 hyperfine transition frequency, and 
7 is the natural width. A continuous wave, low intensity 
diode laser serves as probe of the optical transmission of 
the MOT and also as the excitation source for the time 
resolved experiments reported here. The laser, which has 
a bandwidth ~ 1 MHz, is tuned in a range of several 7 
around this transition. The average light intensity for 
the probe is ^ 1 ^W/cm?. To produce a nearly Gaus¬ 
sian beam profile, the laser output is passed through a 
single-mode optical fiber and expanded to a 1 /e^ width 
^ 8 mm. The probe laser intensity is modulated with an 
acousto optic modulator (AOM), which generates sharp 
rectangular pulses having an on time of 2 /is and an off 
time of 2 ms. The 2 fis excitation pulse is centered in a 90 
/is window during which fluorescence data are recorded. 
The MOT lasers are off during the data acquisition pe¬ 
riod. For the remaining nearly 2 ms, the MOT lasers 
are turned back on to regenerate the ultracold atomic 
sample. An important component of the experiment is 
a synchronized mechanical chopper, which prevents in¬ 
tense fluorescence from the MOT region, present during 
the MOT build up period, from reaching the photomul¬ 
tiplier tube (PMT). The AOM-limited 20 dB response is 
estimated to be less than ^ 60 ns. The probe laser is 
linearly polarized in the vertical direction. 

Fluorescence signals are detected in a direction orthog¬ 
onal to the probe laser propagation and polarization di¬ 
rections. To minimize reduction in the measured light 
polarization, the scattered light is collected in a small 
effective solid angle of about 0.35 mrad, and refocussed 
to match the numerical aperture of a 400 /rm multimode 
fiber. A linear polarization analyzer is used to select the 
detected polarization channel, which we label as parallel 
(II) and perpendicular (T). The fiber output is filtered 
by a 5 nm spectral width interference filter centered at 
780 nm, and passed to a GaAs-cathode photomultiplier 
tube (PMT). The PMT is operated in a photon counting 
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mode, whereby the output is amplified and directed to a 
discriminator and multichannel scalar (MCS) (indicated 
in Fig. 1 as data acquisition and storage). The multi¬ 
channel scaler time sorts and accumulates the data in 5 
ns bins. Timing of the experiment, including the MOT 
and probe lasers, and triggering the start of the MCS 
sweep is controlled by a precision pulse generator slaved 
to a timing signal from the mechanical chopper. 


III. THEORETICAL APPROACH 


We provide here only an outline of our theoretical ap¬ 
proach to description of the time evolution of light scat¬ 
tering in an ultracold atomic gas. Details of this treat¬ 
ment may be found in our earlier papers 0, [lE HOi HH 
The time dependence of the scattered pulse intensity can 
be extracted by using the following calculation steps. Ba¬ 
sically the process is fully described by the scattering 
amplitude defined for each spectral component of the in¬ 
coming and outgoing pulse. Thus, this scattering am¬ 
plitude has to be established as a first step. Such an 
amplitude for the multiple scattering process develop¬ 
ing in an opaque medium includes two important cal¬ 
culation ingredients, namely, the scattering tensor and 
the mesoscopic Green function. The description of these 
theoretical characteristics of the scattering problem for 
the relevant case of an ultracold atomic sample can be 
found in Ref. M]- Roughly speaking the amplitude can 
be constructed for any randomly chosen chain of atomic 
scatterers as a subsequent product of the respective scat¬ 
tering tensors responsible for the deflection of the light 
ray from its coherent forward propagation. The transfor¬ 
mation of the field amplitude inside such a zigzag-path, 
but along the fragments of freely propagation, can be 
described in terms of the Green function formalism. 

In the next step the partial contribution to the time 
profile of the outgoing light pulse is recovered via a re¬ 
verse Fourier transform. This transform should be sub¬ 
sequently evaluated for each order of the multiple scat¬ 
tering. This part of the simulation procedure is the most 
difficult and needs an accumulation of the quite extended 
database for the large number of the contributing spec¬ 
tral components, see [s^ for details. Finally the par¬ 
tial contribution to the instantaneous intensity of the 
light scattered by the selected chain is expressed by the 
squared time dependent partial amplitude. The crucial 
point of the final calculational step comes from the im¬ 
portant physical requirement to the interference contri¬ 
bution associated with the scattering in backward direc¬ 
tion. For this particular scattering channel the partial 
amplitudes responsible for the light propagation along 
two reciprocal scattering path can interfere and the in¬ 
terference survives the configuration averaging. At the 
final step the normal ’’ladder”-type and, in the case of 
the backscattering channel, the ’’crossed”-type amplitude 
products should be averaged with the Monte Garlo tech¬ 
nique for a Gaussian-type distribution of the atomic scat¬ 


terers in the sample. The total outgoing pulse is finally 
expressed by the sum of such partial contributions, and 
the sum is typically converged for the orders of multiple 
scattering ^ bo. 

The described procedure of numerical simulation al¬ 
lows us to include in the numerical routine most of the 
practically important factors such that the hyperfine and 
Zeeman structure of the ground and excited states, the 
spatial inhomogeneity of atomic cloud and its shape vari¬ 
ation. In addition, the dynamical modification of the 
scattering process caused by atomic motion can be traced 
via the time and spectral dependencies of the correlation 
function of the scattered light. Some interesting and im¬ 
portant effects for the light scattering on spin-oriented 
atomic systems were predicted with such a detail numer¬ 
ical analysis, see [2l|. 


IV. RESULTS AND DISCUSSION 

In an earlier paper we reported on the time dependence 
of the linear polarization degree of the multiply scattered 
light. In the measurements, the ultracold atoms are ex¬ 
cited by a nearly rectangular 2 ^s long pulse of linearly 
polarized resonance radiation tuned in the spectral vicin¬ 
ity of the F = 3^F' = 4:^F = 3 transition. The 
measured intensities of linearly polarized scattered light 
in two orthogonal directions (along and perpendicular to 
the direction of the exciting light pulse polarization) may 
be quantified by defining a linear polarization degree as 


Pl 


/||+/U 


(4.1) 


In the formula, /y and I± represent the measured in¬ 
tensities in the lin || lin and lin T lin channels. It was 
shown by Balik, et al. 0 ] that the linear polarization 
degree decays rapidly towards zero on a time scale of less 
than 100 ns. In the present paper, we are concerned with 
the longer time dynamics, for which /y and very nearly 
equal. Our results are then presented solely in terms of 
■^ 11 - 

The time dependence of the scattered light intensity /y, 
for several different detunings from the F = 3 —>■ F' = 4 
transition is shown in Fig. 2. There we see that, after 
a transient growth period of approximately 100 ns, the 
fluorescence intensity reaches a steady state value. Upon 
shut off of the excitation pulse at a time t = 2000 ns, 
the fluorescence decays to a small level on a time scale 
of several hundred ns. As the frequency of excitation 
is shifted away from the resonance transition, both the 
transient build up and decay become more rapid. This 
qualitatively makes sense, for the optical depth of the 
vapor is reduced with detuning from resonance, and the 
radiation thus undergoes fewer scattering events before 
emerging from the optically dense sample. As indicated 
by the solid lines in Fig. 2, the overall time dependence is 
well described by the results of Monte Garlo simulations. 
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FIG. 2: Time-dependence of the scattered light intensity, in 
the /|| polarization channel, for several different detunings 
from atomic resonance. 


FIG. 4: Influence of the laser linewidth on the scattered light 
intensity, in the 7|| polarization channel, for a detuning A = 
-I- 3 MHz. 
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FIG. 3: Semilogarithmic graph of the time-dependence of the 
scattered light intensity, in the /|| polarization channel, for 
several different detunings from atomic resonance. 

The behavior of the intensity after the excitation pulse 
is turned off is seen more clearly in the semilogarithmic 
plots shown in Figure 3. There we see a surprising result: 
the long time decay is nearly exponential and has essen¬ 
tially the same decay rate, regardless of the detuning. We 
point out that, even though this data is for positive de¬ 
tunings of A = 0, 3, 6, and 12 MHz, a very similar effect is 
seen for negative detunings of the same magnitude. It is 
also evident that the Monte Carlo simulations, shown in 
this graph as solid lines, do not account for the effect just 
described. This suggests that additional physics must be 
included to better describe the physical circumstances of 
the experiment. 

One possibility is that the finite bandwidth of about 1 
MHz of the excitation laser results in resonance excita¬ 
tion with the wings of the laser spectral profile. Simula¬ 


tions of this effect for laser widths of 0.01 MHz, 1.0 MHz, 
and 3.0 MHz are shown in Fig. 4. There it is seen that 
for small detunings (A = 3 MHz in this case), there is 
some influence, and the change goes in the proper direc¬ 
tion in comparison with the experimental data. However, 
the bandwidth effect becomes increasingly negligible for 
larger detunings, and so cannot account for the results of 
Fig. 3. 

The most promising possibility is that spectral compo¬ 
nents of the excitation pulse, generated by the relatively 
rapid turn on and turn off of the pulse, sufficiently over¬ 
lap the atomic resonance transition to explain the effect. 
In Fig. 5 we compare the experimental data at a larger 
detuning of A = -|- 12 MHz with simulations for pulses 
having different turn-on and turn-off models and time 
constants. There it is seen, for example, that for an ex¬ 
ponential turn off shape with a time constant on the order 
of a few ns, the data and simulations are in reasonable ac¬ 
cord. As the Fourier components for a rectangular pulse 
drop off with frequency relatively slowly with offset from 
the carrier frequency, this is a suggestive interpretation 
of the observations. We point out that this interpretation 
is supported by the transient turn on behavior, as shown 
in Figure 6. There we see a sharp overshoot, which we 
interpret as beating between the Rayleigh scattered light 
and the Fourier components of the excitation pulse scat¬ 
tered according to the Lorentzian response of the atomic 
resonance. Although the experimental data is too noisy 
to extract the beat frequency, the simulations, on an ex¬ 
panded intensity scale, show that the beat note frequency 
is in fact given by the detuning from resonance. 

Finally, we point out that Labeyrie, et al. have 

observed a similar effect in ®®Rb but in a somewhat differ¬ 
ent fluorescence geometry, sample geometry, and prepa¬ 
ration. In these studies, the Gaussian radius of the MOT 
is 2 - 3 mm, while the peak optical depth bo ranges from 
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FIG. 5: Influence of the temporal shape of the excitation pulse 
on the long time behavior of the scattered light intensity in 
the /|| polarization channel. A = + 12 MHz. 



t (ns) 

FIG. 6: Influence of the temporal shape of the excitation pulse 
on the short-time growth of the scattered light intensity in the 
7|| polarization channel. A = -|- 12 MHz. 

2 to 40. In addition, optical excitation of the sample was 
done with a probe laser smaller than the sample size, in 
contrast to the experiments reported here. Finally, the 
angle of observation of the fluorescence was about 17° 
from the forward direction, in comparison to our mea¬ 
surements where this angle is ~ 90®. They interpret the 
long time constant as being due to redistribution of the 
multiply scattered light by the residual thermal motion 
of the atoms in the ultracold atomic cloud. Reported 


simulations are also in good agreement with their tem¬ 
perature dependent measurements of the decay rate at 
longer times. No data on the turn on portion of the fluo¬ 
rescence measurements were reported. In comparing the 
results of Labeyrie, et al. [^, @ we see that the essen¬ 
tial effect, that the longer time decay is nearly constant 
and independent of detuning, is the same in both ex¬ 
periments. However, the interpretation of the results is 
different. This could well be due to the quite different 
experimental arrangements and conditions for the two 
experiments, under which different physical mechanisms 
may be important. For the experiments reported here, 
it seems that the interpretation based on overlapping 
Fourier components provide a self consistant explanation 
for both the longer time decay and for the shorter time 
turn on transient. However, we should point out that ini¬ 
tially warmer atoms also have a larger overlap with the 
Fourier components of the excitation pulse, and can con¬ 
tribute to the long time tail. This should have a prompt 
effect on the transient buildup and decay, rather than 
an time-accumulating effect such as frequency redistri¬ 
bution, due to thermal motion, in multiple scattering. 
Finally, the data reported in Fig. 5, for example, show 
not ideal agreement with our simulations, even account¬ 
ing for the Fourier components of the excitation pulse. 
This residual difference may be due to thermal redistri¬ 
bution, which has not been taken into account in our 
simulations. 


V. CONCLUSIONS 

In conclusion, we report experimental and theoretical 
results on the long time decay of near resonance radiation 
from an ultracold gas of Rb atoms. The experimental 
results show two features suggesting that Fourier com¬ 
ponents of the excitation pulse play an important role 
in the observed constancy of the scattered light intensity 
at longer times following pulse shut off. Our simulations 
also show that realistic laser spectral widths can influ¬ 
ence the observations for smaller detunings from reso¬ 
nance. Measurements simulations by other researchers, 
on the other hand, have shown that frequency redistribu¬ 
tion due to residual thermal motion of atoms in the cold 
gas plays a significant role in interpreting off resonance 
light diffusion in an ultracold atomic sample. 
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